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TECHNICAL MEMORANDUM 


RADIATIVE AND FREE CONVECTIVE HEAT TRANSFER 
FROM A CONTAINERLESS SPHERE 

INTRODUCTION 


Containerless processing of copper in an electromagnetic levitation 
device uses the natural heating of the sample (generated by the induced 
eddy currents) and cooling by a helium-argon gas mixture. Such a 
system has been in operation in the Space Sciences Laboratory, Marshall 
Space Flight Center, and significant undercooling during solidification 
has been noted. Temperatures on the order of 2000 K were known to be 
reached when approximately 20 torr of argon were present in the levita- 
tion chamber , and small amounts of helium (approximately 100 torr 
pressure) were known to reduce the copper sphere's temperature by 
more than 1000 K. As the coolant gases were admitted to the evacuated 
chamber through a valve at a large relative distance to the sphere, heat 
transfer could be considered to be accomplished by radiation and free 
convection . 

In support of this experimentation, a relatively simple mathematical 
model for the prediction of heat transfer from the sphere was derived. 
Primary conditions to be considered were the nature of heat loss, the 
high temperatures and low pressures involved, and the restrictions on 
direct experimental measurements of various parameters (sphere diameter, 
emissivity and temperature, and gas temperature and pressure were the 
only measurable quantities). Parameters to be varied in the calculation 
of temperature were power input to the levitation coil, emissivity and 
diameter of the sphere, and gas pressures. 


DETERMINATION OF HEAT TRANSFER FORMULA 


Power input (PI) to the sphere is equal to total heat loss through 
radiation and free convection: 


PI = Q , + Q 

^rad ^conv 


The heat transfer due to radiation (Q rad ) is 
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where 

■’f - Stefan-Boltzman constant ^ Q.56G97 x l(f 8 W/m 2 -K 4 
s “ cmissivity of sphere 

o 

A “ surface area of sphere, cm 
- sphere temporuture, K 
T ~ ambient gas temperature, K , 

The heat transfer clue to free convection (Q conv ) is 
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where h ~ heat transfer ratio for free convection [ 1 ] or, putting h in 
terms of the Nusselt number, Nu, 

Q o = A(T 1 - T q ) J 5 Nu 


where k - thermal conductivity of coolant gas and D = sphere diameter. 

Yuge [2] empirically determined the Nusselt number for a sphere 
in air with heat loss due to free convection to be 

Nu = 2 + 0.39 Gr 1/4 for 1 < Gr < 10 5 


where Gr = Grashof number. For Grashof numbei’s much less than one, 
Yuge recommended Mahony's equation: 

Nu = 2 + 0 (Gr 1/2 ) 


where 0 = order of magnitude. By fitting a line to the lower portion of 
Yuge's data, 0 was found to equal 0.25. In his work, Mahony dealt 
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mainly with Grashof numbers between 10~ 9 and 1(T 9 [3], Because the 

- 3 

Grnshol numbers encountered fell between 10 and one, both formulas 
wore tried. No significant differences in results were found. Yuge's 
equation was used because of the greater documentation and certainty 
in the Grashof coefficient. 

The Grashof number was calculated in the same manner as Yuge: 


Gr = gD 3 AT 

T Vp 
of 

where 

g = acceleration due to gravity 
D = sphere diameter 

AT = sphere temperature T^ - gas temperature T q 

T q = gas temperature; 1/T = heat transfer coefficient (3) 

T 1 + T o 

Vj, = kinematic viscosity at film temperature g 



where 

= gas density at film temperature 
Uj, = dynamic viscosity at film temperature. 


PROPERTIES OF COOLANT GASES 


Because no direct measurement of gas properties (with the excep- 
tions of pressure and temperature) was possible, methods for estimating 
the thermal conductivity, dynamic viscosity, and density were either 
derived or adapted. All properties were found to be highly temperature 
dependent, and the temperature difference T^ - T between the sphere 
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mid the surrounding’ gas was very largo. Therefore, all properties 
wore evaluated ut the film temperature T f - (T. - T )/2, ns 
recommended by Yuge. 1 to 

In the. low pressure range found in the experiments, thermal con- 
ductivity k is not pressure dependent [4]. A temperature dependency 
3 /4 

of Tp was determined empirically from data given in McAdams [1]: 


/ rp \ 0.75 

k = k Q (0.0173) ( 2 Y 3 ) W/cm K 

whore 

k = thermal conductivity of specific gas at 273 K in Btu/hr *ft~°F 
0.0173 = conversion factor from Btu/hr-ft-°F to W/cm K. 


This formula wus weighted to provide more precise correlation with 
McAdams' data for higher temperatures and helium, the major component 
gas . 


Since mixtures of helium and argon do not involve rotational and 
vibrational degrees of freedom, Brokuw's simple empirical method for 
calculation of thermal conductivity [5] could be used. Brokaw utilized 
the fact that thermal conductivities for mixtures fall between the values 
found by simple mixing of the component conductivities and inverse 
mixing of the components: 
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where x^>X 2 " m °l e fractions of component gases. Brokaw represented 
the mixture conductivity as a combination of simple and inverse mixing: 


k m ~ qk SM + (1 “ q)k RM 


where 


k SM y l k l + X 2 k 2 
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and 



Brokaw then fit tills curve to experimental data of numerous binary gas 
mixtures, including helium-argon, for various molar fractions to deter- 
mine the value of q, The value of q was found to vary with the ratios 
of the gases [5]: 

mole fraction 0 .1 .2 .3 ,4 .5 .6 .7 .8 .9 .95 1.0 

q .32 .34 .37 .39 .42 .46 .50 ,55 .61 .69 .74 .8 

Brokaw' s method was found to have approximately the same 
accuracy as other methods studied and to be much sampler in form. 
Therefore, Brokaw's method for gas mixtures was used in combination 
with the author's formula for individual gases to estimate thermal con- 
ductivities of the coolant gases. 

Because the dynamic viscosities and the densities of the gases 
were needed in the calculation of the Grashof number, methods to calcu- 
late these properties were also needed. For the component gases, 
dynamic viscosity g was taken to be 


rtf 

p = '‘o-llvn Poiso 


where y = dynamic viscosity of specific gas at 373 K (P) , again by 

fitting a curve to McAdams' data [1]. A pressure dependency was not 
incorporated due to the high temperature and low pressure values to be 
used. 


Brokaw's method for calculation of viscosity of gas mixtures was 
used, as recommended by Reid, Prausnitz, and Sherwood. Brokaw's 
method is based on Sutherland's approximation 
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where 

Sy s Sutherland constant 

A., - a complicated function of molecular weights of the component 
1J gases . 

A., can be taken from a chart [6] or a graph [4] and, for mixtures of 
nonpolar gases, Sj. can be set equal to one, Thus, for mixtures con- 
taining nonpolar gases, 
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This mixture formula has its best results for mixtures of inert gases; 
overall ex’rors for nonpolar gas mixtures ranged from 0.6 percent to 
2.5 percent error [6], Brokaw’s method can also be applied at any 
temperature, a vital feature for this project. 

The densities of the gases were calculated by the formula 


P 



where Tj, is in K, P is in torr, and M = molecular weight of gas. The 
method of simple mixing was used to calculate the mixture density: 


m 


= X 


l p l 


+ X 2 P 2 
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THE COMPUTER PROGRAM 


A FORTRAN program was designed to simplify calculations and 
extend the range of applicability to experimentation. Appendix A pre- 
sents the program nomenclature, including a cross-reference to the 
nomenclature used in the text of this report; Appendix B presents a 
program listing. 

In the program, nested DO loops are used to allow for variances 
of sphere emissivity EP, sphere diameter D, and power input PIN (Set 
Experimental Parameters section) in addition to argon pressure PI. Two 
separate DO loops are provided to change gas pressure (Inareaoe 
Pressure and Decrease Pressure sections) , simulating the experimental 
input and pump down of gas. The Increase Pressure loop uses the 
incremented mole fraction to calculate helium pressure P2 for constant 
PI while the Decrease Pressure loop steadily decreases the total pressure 
PRTTL, holding the molar ratio of the gas mixture constant. Each 
pi’essure loop includes an implied loop within which the gas properties 
and sphere temperature are calculated (iterate to Satisfy Temperature 
Equation sections). The gas properties arc determined by solving the 
equations presented earlier in this report; the temperature is found 
using an iterative technique wherein the basic heat transfer is solved 
for sphere temperature T by using an estimated sphere temperature . 
T^ is initialized at RADT, the sphere temperature should only radia- 1 

tive heat transfer be present. The T^'s are revised as the implied loop 

is repeated until the difference |T - TjJ is less than an input tolerance 

(in this case, 10 K). An iteration check (FLAG, FLAGM) is also pi’o- 
vided to check for bad convergence. 


RESULTS AND CONCLUSIONS 


As pressure is increased, convection takes on an increasingly 
important role in cooling the sphere. Once the helium enters, convec- 
tion becomes dominant. Varying the amount of argon initially present 
does not significantly change the sphere temperature (Figs. 1 and 2). 

The percentage of helium present in the mixture, and not the total 
amount of helium, determines the sphere temperature: 90 torr of helium 

added to 10 torr of argon changes the sphere temperature from 1577 K 
to 637 K, while 360 torr helium added to 40 torr argon effects a change 
from 1519 K to 618 K. This phenomenon is due to the fact that the 
Grashof number is the only pressure- dependent term in the heat transfer 
equation. As expected, emissivity produced the majority of its effect at 


low pressure, whore Hod T is proportions to omissivity (Fig. 3). 
Sphere temperature is directly proportional to power input (Fig-, 4) and 
inversely proportional to sphere diameter (Fig. 5) . Once convective 
cooling dominates, an increase of 1.5 W in power is comparable to a 
decrease of 0.135 cm in diameter. During radiation's dominance, the 
power input change will not be as effective as the change in diameter. 


Small fluctuations in the value of k , will effect great changes in 

nux 

the values of Comb T, while Bad T will remain constant. It is believod 
that the major source of error in this program is in the value of k j x , 

As temperatures increase, the estimation for individual thermal conduc- 
tivities will tend to bo high, This deviation will decrease the apparent 
pressure effect on sphere temperature Comb T, 


Correlation for tho Ineveaoe Premiere portion of the program with 
experimental results was found to bo good: tho variance between sphere 
temperature at 20 torr argon and 200 torr total was 2,7 percent, Tho 
actual temperature values wore approximately 5 percent higher than 
expected: however, tho error in calculation of thermal conductivity would 
predict a higher power input (and thus temperature) for the same cooling 
effect on the sphere. 


Results for the Dcereaac Premiere portion of the program were 
discouraging. The increase in temperature is only approximately one- 
third that found experimentally, This is believed to be due to the 
dominance of the conductivity term, which depends upon the ratio of the 
gases rather than the total amount of gus present (Fig. 6). The small 
temperature change shown illustrates the temperature dependence on the 
Grashof number. 
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NOTE TO FIGURES: 


Because total heat lass did not vary significantly with the 
total amounts of gas present but rather with the percentage of helium 
present (Figs, 1,2), Figures 2 through B are plotted as temperature 
versus percentage helium. AH lata (except where noted) wore taken 
from rims with 20 torr argon initial pressure. The following is an 
approximation of results for runs beginning with 10 or 40 torr argon: 




Total Pressure 


He, Percent 

To™Foin~Ar"" 

20 tonTAr 

40 torr A? 

0 

10.0 

20.0 

40.0 

10 

11.1 

22.2 

44.4 

20 

12.5 

25.0 

50.0 

30 

14.3 

28.6 

57.1 

40 

16.7 

33.3 

66.7 

50 

20.0 

40.0 

80.0 

60 

25.0 

50.0 

100.0 

70 

33.3 

60.7 

133.3 

80 

50.0 

100.0 

200.0 

90 

100.0 

200.0 

400.0 
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© EMISSIVITY = 0.10 
a EMISSIVITY - 0.28 



POWER INPUT 

2.5 WATTS 


SPHERE DIAMETER 

0.385 CM 


ARGON PRESSURE 

20 torr 


INTRODUCTION OF 
20 torr ARGON 


PERCENTAGE HELIUM (%) 


Figure 3. Sample temperature as a function of percentage helium 

for various emissivities. 
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300 A 


EMISSIVITY 
POWER INPUT 
ARGON PRESSURE 


= 0.16 
= 2.5 WATTS 

= 20.0 tprr 


200 A 


I { 20 40 60 80 100 

INTRODUCTION PERCENTAGE HELIUM (%) 

OF 2Qt ARGON 


Figure 5* Sample temperature as a function of percentage helium 
for various sphere diameters. 


THERMAL CONDUCTIVITY OF GAS MIXTURE 1C MIX X10* 3 (POISE) 





NOTE: BOTH THERMAL 
CONDUCTIVITY AND GRASHOF 
NUMBER ARE FUNCTIONS OF 
PRESSURE (EITHER TOTAL, OR 
AS REPRESENTED BY PERCENTAGE 
HELIUM) AND TEMPERATURE. 
THUS, THE VARIANCES SHOWN 
HERE ARE NOT DUE SOLELY TO 
THE VARIANCE OF PERCENTAGE 
HELIUM. INCREASE IN 
PERCENTAGE HELIUM 
IMPLIES A DECREASE IN 
TEMPERATURE. THIS 
TEMPERATURE DEPENDENCY 
IS INCORPORATED IN THE 
CURVES SHOWN. 


Figure 6. Thermal conductivity and Grashof numbers as functions 

of percentage helium. 
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APPENDIX A 


NOMENCLATURE 


Report 

Program 

Definition 

A 12' A 21 

A12,A2l 

Dimensionless parameter used in calculation of 
gas mixture viscosities, function of the molecu- 

/ M i w 2 \ 

lar weight ratios ( j 

" 

ARGT 

The argument of T used in calculation of inter- 
mediate values of sphere temperatures T - 

ARGT 1/4 (K 4 ) 

~ 

AVDT 

Absolute value of intermediate sphere tempera- 
ture value difference ]T - T.J (K) 

Comb T 

COMBT 

Sphere temperature due to combined radiative 
and free convective heat transfer (K) 

D 

D 

Diameter of sphere (cm) 

5 

EP 

Emissivity of sphere 

FLAG 

FLAG 

Iteration counter for temperature calculation 

FLAGM 

FLAGiVI 

Maximum number of iterations allowed in calcu- 
lation of sphere temperature 

Gr 

GR 

Grashof number for gas mixture, dimensionless 
ratio of buoyant to inertial forces 

k 

o 

K10.K20 

Thermal conductivity of gas (1,2) at 273 K 
(Btu-hr-ft-°F) , used in calculation of thermal 
conductivity of gas (W/cm-K) at film temperature 

k l ,k 2 

I< 1 , K 2 

Thermal conductivity of gas (1,2) (W/cm-K) 

k m 

KMIX 

Thermal conductivity of gas mixture (W/cm-K) 

M 

Ml, M2 

Molecular weight of gas (1,2) 

y 

*0 

MU 10 ,MU20 

Dynamic viscosity of gas (1,2) at 373 K (P) 
used in calculation of dynamic viscosity of gas 
(1,2) at film temperature (P) 
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NOMENCLATURE (Concluded) 


p^r^yi r ** " 


Report 

Program 

Definition 

M 

MU1,MU2 

Dynamic viscosity of gas (1,2) at film tem- 
perature (P) 

p m 

MUM IX 

Dynamic viscosity of gas mixture at film tem- 
perature (P) 

P 1’ P 2 

P1,P2 

Amount of pressure of gas (1,2) (torr) 

~ 

PCNTl,PCNT2 

Percentage of gas (1,2) in mixture (percent) 

PI 

PIN 

Power input to levitation coil (W) 

P ttl 

PTTL 

Total gas pressure during increase of pressure 
(torr) 

PRTTL 

PRTTL 

Total gas pressure during pump down (torr) 

q 

Q(D 

Dimensionless parameter used in calculation of 
thermal conductivity of gas mixture, function 
of molar fraction of light gas 

Rad T 

RADT 

Sphere temperature when only radiative cooling 
effects are considered (K) 

p l’ p 2 

RM01,RH02 

Density of gas (1,2) 

p m 

RHOMX 

Density of gas mixture 

T 

T 

Intermediate sphere temperature value used in 
determination of sphere temperature (K) 

T o 

TO 

Ambient gas temperature (K) 

T 1 

T1 

Intermediate sphere temperature value used in 
determination of sphere temperature (K) 

T 

X f 

— 

T 1 + T o 

Film temperature ^ 

— 

TOL 

Tolerance between intermediate and final sphere 
temperature values (K) 

Xi*X 2 

XI, X2 

Mole fraction of gas (1,2) 
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•090 RH02 * ( 3 ,21E-S*M2tP2)/(Tl+T0 ) 

•• ■'•■ RHOMX * X 1 * RHO 1 + X2 * RHQ2 

>■>>*•;; MUt ■ 3 66£-2 + MUlO + <TmO>**,5 

•■■•'?» 3 HU2 * 3 6GE-2*MU20*<T1 + TQ>»+ 5 

>"/94 MUMIX » <<XUMUl V<Xl+A12*<MUl/HU2;i**.3»X2>> + <<X2*Mll2V<X2+A21* 

'• 1 < NU2/MU1 )* * . 5 *X 1 ) > 

•■ 09 ,. K 1 * 1 53E-4*K10t< < T14T0)M 75) 

■••'-'* K 2 * 1 53E-4*K20»<<Tl+TQ)** 75' 

ov'*« KM1X a CU 1>*<X1*K1+X2*K2H<1 ,»(l(I))*(l . A < X 1 XK 1 > + < X 2/K 2 > ) > 

■■•‘• l IF MI2- 3) 65, 65,93 

i«'<H OS ARGTol <PIN-<6 28*KHIX*D*<Tl~T0>)-< < 6 . 86 *K H l X *D ♦ ♦ 1 . 7S»RH0MX*». 5* 
•■loj 1 >T1-T0)*M 25>X(T0t* ,25*HUMIX** . 5) ) V < 1 78 E- 1 2* EP* D* »2 ) > 

•'*,'• ,' t + TO » * 4 

•!0’i I F 1 A R G T ' 70,70,90 

■<*.04 ;■; IF, FLAGH-FLAG ) 1 40,14 0, 80 

-i.'*-. hv T 1 * T 1 - 100 

• * 0 1 < G 0 T 0 6 v 

4- T n ARGT * * 25 
GO TO H 

■ : < 17 T °« < P I IV * 1 7SE-12>EP-+P**2+<TtM4-T0M4 ))~U 6 8 6* K M I X ♦ 0 ♦* 1 . 7 5* 

■*•11 ■ ■ F H 0 N it <■ *■ ,5*< T 1 » TO > ♦ * I 2 S >A T 0 * » 2 5 * MUM I X ♦ > . 5 > ) > / ( 6 28*KMIX*P)> 

■■; ! , + T 0 

• 1 ' ? 4f. A VP T « CABSf T-Tl ) 

V* J i IF** TOL-AVDT 1 1 00, 1 60, 160 

>, t ! 4 too IF FLAGM-FLAG' 1 40,1 40, 1 10 

-11 5 I 1 <* 11 n < T+Tl 

• 5 • 1 F* T 1 - TO * 1 20, 120, 130 

•> ! f ;,‘:0 Tt « TO + too 


•If "(I GO to 60 


00 320 1 EP* 10 * 26 ■ b 
EP « FLOAT. 1EP ,,*100 
00 310 IP*250.3B5,135 
0 * FLOAT < I 0 '<*1000 
00 300 I P 1 N * 1 0 40, 1 5 
PIN * FLOAT! 1 P IN VI 0 
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j I H 

"U’V 
*52 1 
01 £2 
•H23 
0 124 
0125 

0 12 W 
012? 
OUQ 

0129 

0130 

0131 
«>132 

0133 

0134 

0135 
013$ 

0132 

0138 

0139 

0140 
014 1 

0142 

0143 
0 t 4 4 
0 !45 

0146 

0147 
014U 
0 149 
ul 30 
M 1 5 1 

0152 

0153 

0154 

0155 

0156 
0 157 

0158 

0159 

0160 
0161 
0162 

0163 

0164 

0165 

0166 
016? 
0168 

0169 

0170 

0171 
0 172 

0173 

0174 

0175 

0176 
0 177 
0 J 78 


C CALCULATE AND WRITE FI HAL DATA 

i* 

140 WRITE v 6 * 150 > T, T l , RADT 

150 FORMAT v ,'* BAD CONVERCEHCE THE LAST TWO VALUES OF T WERE ' • F6 0, 
i * AND *,F6 0,' RAOT ■ '.FM7) 

GO TO 180 

160 PTTL ■ PI f P2 
PCHT2 > X2* 100 . 

CR * C 960 ♦ D**3*RHQNX**2*!Tl«T0>>,'<T0*MUNIX**2> 

WRITE! 6. 170) EP , P I H , 0 » P T TL , P l , P2, P CH T2 , RAPT , T . OR , KM IX 
170 F0RNAT!F5.2,F11. liFll 3,Ft31,3Ft2.1,F13.0,F12.0,El33,El3.3> 

180 CONTINUE 

C 

C 

C DECREASE PRESSURE 

C 

C 

URITE <6, 185)^ 

7/* DECREASE PRESSURE* 7/3X, *EP * , 8X , "P IN * , 9X , * D * , t OX, 
I'PTTL*,8X.*PRES1*,?X. 'PRESS’ . 7X . » PCH T2 * . 7X. " RADT * . 8X. * CONST*. 10X. 

< ‘GR* . tOX, *KNI X "7 ) 

C 

C SET CONSTANTS AND INITIALIZE VARIABLES 
C 

XI ■ 1 

X 2 ■ 9 

PCNT2 ■ X2 * 100 
PRITL ■ PTTL 
T 1 ■ RADT 
DO 290 H • 1 . 13 
PRES 1 • X 1 * PRT TL 
PRES2 ■ X2*PRTTL 
FLAC « 0 
C 

C ITERATE TO SATISFY TEMPERATURE EQUATION 
C 

190 FLAC «• FLAG * 1. 

RHOl » < 3.21E-S*Hl*PRESl >/<Tl+TO> 

RH02 ■* <3.21E-5*N2*PRES‘2)7!T1+T0> 

RHOMX ■ X1*RH01+X2*RH02 

MU1 •> 3>66E“2*MUI0*<T 1 + TO > *♦ . 5 

II U 2 * 3. 66E-2*MU20*<Tl+T0>** .5 

MUM IX * !<Xl*MUl)7!Xl+A12*<HUl7HU2)**,5*X2)>K<X2*HU2>7!X2'f'A2l* 

I < MU27MUI )** ,5*Xl )> 

K 1 ■ 1 53E-4*K10*<!mTO>**,75> 

K2 « 1 .53E-1*K20*! < T 1 + T 0 )* * , 75> 

KMI X ■ Q< I>*<X1*K1+X2*K2> + <1 , -0< I >)*<1 . 7! ! X 1 7K 1 >K X27K2> ) > 

T » << PIN-< 1,?BE-12*EP*D**2*<T1**4-T0**4>)‘-<!6,86*KNIX*D**1 .73* 

! RHOMX** 5*<T1”TO>**I.25>7<TO*»,23*HUMIX**,5>>)7<6,28*KN1X*D>) 

I +1 0 

A VP 1 r PARS! T- T 1 ) 

IK 1 01 -AVD) ) 200,260.260 
200 IF' 1 FLAGH-Fl AG ) 240,240, 210 
210 11 » <H1) >72. 

IF! T 1 - 1 0 > 220,220,230 
220 T1 *> TO ♦ 100. 

230 CO 10 190 

C 
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REPRODUCIBILITY OF THE 
ORIGINAL PAGE IS POOR 


OlPO 
0161 
0182 
0181 
0184 
0 185 
0160 
018',* 
0188 
0168 
0190 
0 191 
019? 

0193 

0194 

0195 
0 196 
019? 

0196 

0199 

0200 
020 1 
0202 
»>D03 
•>?0<t 


» uacumr ahp hr me fimai pata 

e 

240 8RIIM6.250) T.t» 

2£<0 FORMAT (/* BAP CONVERGENCE 1 « »,(■*. 8, • Tl « *,F6 0/J 
CO 10 260 

260 fR * < 960 «PM3*RH0HXM2*ni-I0))/<T0*MUHlXo2> 

V0MC(«.27O) f P»PlK<D<FRTH.PRESI»PRE92»PCNT2»ftAPT»T>CR»KHIX 
25 0 FOiVH AT (F 5 2, F 1 1 . 1 , F 11 - 3. F 1 3 l,3Fl2.l,Fi3 0,F12.0,E15 3,613.3) 

r 

C CHECK IF PUN? P08N IS COHPUTEP 

C 

PITH » PRTU-30. 

1 F< PRT TL LE.O, ) CO 
280 CONTINUE 
265 URITEU.P90) 

290 fQPhrtT<mX,'***«** 

!**♦**»****» ♦**?/> 

300 CONTINUE 
310 CONTINUE 
320 CONTINUE 

URJTE< 6.330 ) 

330 FORMAT < I Hi ) 

340 CONTINUE 
STOP 
EHO 
EHPI 


TO 205 
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APPROVAL 

RADIATIVE AND FREE CONVECTIVE HEAT TRANSFER 
FROM A CONTAINERLESS SPHERE 

By Karen Johnson 


The information in this report has boon reviewed for technical 
content, Review of any information concerning Department of Defense 
or nuclear energy activities or programs has been made by the MSEC 
Security Classification Officer This report, in its entirety, has been 
determined to bo unclassified. 
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